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ABSTRACT 

In view of the assumption that any planetary system is likely to be composed of 
more than one planet, and the multiple planet system with a large mass planet has 
more chance of detailed follow-up observations, the multiple planet system may be 
an efficient way to search for sub- Jovian planets. We study the central region of the 
magnification pattern for the triple lens system composed of a star, a Jovian mass 
planet, and a low- mass planet to answer the question if the low- mass planet can be 
detected in high magnification events. We compare the magnification pattern of the 
triple lens system with that of a best-fitted binary system composed of a star and a 
Jovian mass planet, and check the probability in detecting the low-mass secondary 
planet whose signature will be superposed on that of the primary Jovian mass planet. 
Detection probabilities of the low-mass planet in the triple lens system are quite similar 
to the probability in detecting such a low-mass planet in a binary system with a star 
and only a low- mass planet, which shows that the signature of a low- mass planet can be 
effectively detected even when it is concurrent with the signature of the more massive 
planet, implying that the binary superposition approximation works over a relatively 
broad range of planet mass ratio and separations, and the inaccuracies thereof do 
not significantly affect the detection probability of the lower mass secondary planet. 
Since the signature of the Jovian mass planet will be larger and lasts longer, thereby 
warranting more intensive follow-up observations, the actual detection rate of the low- 
mass planet in a triple system with a Jovian mass can be significantly higher than 
that in a binary system with a low-mass planet only. We conclude that it may be 
worthwhile to develop an efficient algorithm to search for 'super-Earth' planets in the 
paradigm of the triple lens model for high-magnification microlensing events. 

Key words: gravitational lensing — planetary systems — planets and satellites: 
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1 INTRODUCTION 



To detect and characterize extrasolar planets various tech- 
niques have been employed, includ ing the radial veloc- 
ity technique dMavor fc Quelo j [i995h . the transit method 
JCharbonneau et al. 20071 ). direc t imaging (jChauvin et al.l 



120041 ). pulsar timi ng analysis ([W olszczan & Frail 1992), 



and microlensing ( Bond et al.l |2004|; lUdalski et al.l [200 
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ISumi et al.ll2010l ). Compared to other techniques, the mi- 
crolensing method has an important advantage of being ap- 
plicable to planets to which other methods are generally in- 
sensitive; the microlensing technique is sensitive to detecting 
low- mass planets and cool planets, or even free-flo ating plan- 
ets (|Bennett fc Rhiell2002l : lHan et al.ir2004l2005l ). This sen- 
sitivity is so important for testing the core accretion theory 
of planet formation which predicts that the dominant plan- 
ets in any planetary system should form in the vicinity of the 
' snow line', which is located at a few AU f rom the host star 
llLaughlin. Bodcnhci mer. fc Adams! I200I : llda. fc Linl [20051 : 
iKennedv. Kenvon. fc Bromlevll2006l ). 



When a microlensing event occurs by a planetary 
system, the planetary signal appears as a short-duration 
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perturbatio n to the standard light curve induced b y the 
lens star dMao fc Paczvnskil Il99ll ; iGould fc Loebl Il992l ; 
iBennett fc Rhielll996l ). The planetary lensing signal induced 
by a planet with a mass of Jupiter lasts for a duration of 
~ 1 day and of ~ 1.5 hours for a mass of Earth. There- 
fore, the discovery of a terrestrial planet would only be 
possible by high-cadence anomaly monitoring. In fact, the 
detection of a significant number of terrestrial extra-solar 
planets requires well-coordinated efforts involving a net- 
work (e.g., MOA-II: Sumi et al. 2010; OGLE-IV: Udalski 
et al. 2005). MOA-II has reported two low- mass planets 
with their survey, MOA-2007-BLG-192Lb (Bennett et al. 
2008) and OGLE-2007-BLG-368Lb (Sumi et al. 2010), and 
is preparing one, MOA-2009-BLG-266Lb (in preparation). 
The planetary perturbation occurs when the source star 
crosses the caustic or passes close to it. Caustic-crossing 
events cause conspicuous double-peaks over the smooth light 
curve induced by a lensing star. However, the perturba- 
tions due to caustic-crossing occur without any prior warn- 
ing so that current microlensing follow-up observations are 
focusing on high magnification events for the sake of practi- 
cality. For high-magnification events, the source trajecto- 
ries always pass close to the perturbation region around 
the central caustic induced by th e planet and thus timing 
can be predicted fairly accurately dGriest fc SafizadehJ | l998 ; 
Han fc Kimll200ll;lBond et alj|2002t iRattenburv et alJbool; 



Yoo et al.ll2004D . lGriest fc Safizadehl ||1998I ) have shown that 
planets with masses as low as IOMq could be detected with 
significant probability in events with magnifications ~ 50 by 
monitoring the peaks of the events with a photometric preci- 
sion of ~ 1 per cent. Rhie et al.(2000) have first showed that 
high magnification events are sensitive to low-mass (Earth- 
mass) planets (For the lower limit of the planet mass most 
recently reported, see Yee et al. 2009). 

Since the discovery of the first extra sola r planet or- 
biting a solar-type star |Mavor fc Qucloz 19951). the Extra- 
solar Planet Encyclopedi^B lists 429 entries, including 45 
multiple planetary systems, as of the 09th of February in 
2010. The detectable mass of exo-planets is extending be- 
low the IOMq regime, wi th the discoveries o f Gliese 876d 
at a mass of ~ 7.5M jRivera et al.ll2005l ). three plan- 
ets aroun d HP 40307 with m asses of ~ 4.2, ~ 6.9 and 
~ 9.2M0 llMavor et alj|2009al). an d Gliese 581e with a mass 
of ~ 1.9A% jMavor et al.ll2009bh using the radial-velocity 
technique, as well as OGLE 2005-BL G-390Lb detected by 
microlensing at a m ass of ~ 5.5M (|Beaulieu et al.ll2006l : 
IBennett et alj 12008V Many of the discovered 'super-Earth' 
planets have been revealed through the close re-examination 
of planetary signals that already proved the existence of 
their big brother. Considering that a planetary system is 
likely to be composed of more than one planet, this kind of 
strategy to find low-mass planets in the multiple planet sys- 
tems may become an efficient way to search for terrestrial 
planets in the sense that it is easier to detect subtle signals 
when one knows where to look for. 

The magnification pattern due to the triple-lens systems 
is known to be well approximated by the superposition of 
the magnifications due to individual planets for planetary 
caustics (Han et al. 2001) and for central caustics (Han 
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2005). Therefore, we may expect that the detection prob- 
ability of low-mass planets in a triple system with a Jovian 
mass planet plus a low-mass planet should be quite similar 
to the detection probability in a binary system with only 
a low-mass planet. However, this is known only for a lim- 
ited range of planet mass ratio and separations (Han 2005), 
and also the effect of the deviation from the superposition 
approximation to the detection probability is not known. 
Moreover, it is not clear if the smaller deviation in the mag- 
nification due to a low-mass planet superposed on the larger 
deviation due to a more massive Jovian mass planet can be 
effectively detected in high magnification events when the 
'adjusted' binary lens model is fitted to incorporate the ad- 
ditional deviation from the low-mass planet, possibly remov- 
ing the signature of the low-mass planet. So in this paper, 
we study the triple lens system (a lens star, a Jovian mass 
planet, and a low-mass planet) for a broader range of plan- 
etary masses and separations than in previous studies, mo- 
tivated by the fact that all planets in the lensing zone will 
substantially affect the central caustics and thus the exis- 
tence of the multiple planets can be inferred by analyzing 
additionally deformed anomalies i n the light curve of hig h 
magnification microlensing events |Gaudi et al.| [l998. 2008). 
We further calculate the probability to detect the low-mass 
planet in the high-magnification events of a triple lens sys- 
tem and compare that to that of a binary lens system (a lens 
star and a low-mass planet). Since the dependence of the size 
of t he central caustic on the planet /star mass ratio q is lin - 
ear (|Griest fc Safizadehl 19981 : IChung et ail2005l : lHanll2006T ). 
detecting signals of Earth-mass planets with q ~ 10 -5 from 
the planet-search strategy of monitoring high-magnification 
event of the binary lens system is observationally challengin g 
(e.g., IBennett fc Rhi3ll996l . 120021 : iGaudi. fc Sacked I2000T ). 
Therefore, it should be crucial and timely to address the 
question if the detection probability of a low-mass planet 
for high-magnification events of the triple lens system is 
comparable to that of the binary lens system. Answers to 
this question may well have implications on interpretations 
of the light curve of high-magnification microlensing events. 
Indeed, recently there have been efforts to re-analyze the 
observationa l data along the line of this point. For example, 
iKubas et al.l |2008l ') have attempted to estimate the detec- 
tion probability for the secondary planet by re-analyzing the 
observational data of OGLE 2005-BLG-390. Their assump- 
tions may be justified for their purposes in that the impact 
parameter of OGLE 2005-BLG-390 is large, i.e., u = 0.359. 
On the contrary, in this paper we specifically concentrate on 
the high magnification microlensing events (uo S 0.01). 

In §2, we begin with a brief description of the multiple 
lens systems. In §3, we construct the fractional deviation of 
the magnification map induced by the secondary planet. In 
§4, we present the detection probability and compare the 
triple lens case with the binary lens case. In §5, we conclude 
with a summary of our results. 



2 MULTIPLE LENS SYSTEMS 

When a source star is gravitationally lensed by a point-mass 
lens, the source is split into two images with the total mag- 
nification given as a simple analytic form of 



DETECTION PROBABILITY OF A LOW-MASS PLANET FOR TRIPLE LENS EVENTS 3 



A = 



vf ± 2 



(1) 



u\/u 2 + 4 

where u is the lens-source se paration normal ized by the an- 
gular Einstein ring radius (|Paczvhskil Il986h . The angular 
Einstein ring radius is related to the physical parameters of 
the lens by 



Oe 



4GM 



\D ol Dos J ' 



(2) 



where M is the mass of the lens and D \ and D os repre- 
sent the distances from the observer to the lens and source, 
respectively. 

If the lens star accommodates a planetary companion, 
the latter may further perturb the image and so changes 
the magnification. When a lensing event is caused by the 
multiple lens system, locations of the individual image are 
obtained by solving the lens equation expressed in complex 
notations by 



( = z + 



N-l 

E — 

3=0 



(3) 



where qj 's and Zj 's represent the mass fractions of individual 
lenses (j = for the central lens star, j = 1 for the primary 
planet, and j — 2 for the secondary planet and so on) such 
that , qj = 1 and the positions of the lenses, respectively, 
£ = f + irj and z — x + iy are the positions of the sourc e and 
imag es, and z denotes the complex conjugate of z (Witt 
ll990T ). Note that all these lengths are normalized by the 
combined Einstein ring radius, which is equivalent to the 
Einstein ring radius of the single lens with a mass equal to 
the total mass of the system. The total magnification is the 
sum of magnifications of the individual images, A — Y^kAk- 
The magnification of each image evaluated at the position 
of each image Zk, is given by 

(4) 



.4 



|det J\ Zk 



where 



det J : 



(•») 



dz dz 

The fundamental difference in the geometry of the mul- 
tiple lens system from that of the single point-mass lens 
system is the formation of caustics. The caustic refers to the 
source position on which the magnification of a point source 
event becomes infinity, i.e., det J — 0. The set of caustics 
forms closed curves. For example, for the planetary lensing 
system composed of a single star and a single planet, there 
exist two sets of disconnected caustics: One 'central' caus- 
tic is located close to the lens star, while the other 'plan- 
etary' caustic(s) is (are) located away from the lens star. 
The number of the planetary caustics is one or two depend- 
ing on whether the planet lies outside or inside the Einstein 
ring. When the planetary companion is close to the Einstein 
ring, the planetary and central caustics merge into a single 
'resonant caustic'. The central caustic plays a crucial role 
in current microlensing planet searches, as stated earlier, in 
that the central caustic due to the planet has a high prob- 
ability of perturbing the light curve of high-magnification 
event, therefore detectable as the planetary signal. For the 
triple lens case, the lens equation becomes a rather compli- 
cated tenth order polynomial equation that has 4, 6, 8, or 



10 solutions, corresponding to physical images depending on 
the configuration of the lens system and the location of the 
source (Rhie 2002) . In the triple lens system, patterns of the 
planetary caustics may well be described by the superposi- 
tion of th ose of the single-planet systems (|Han et al.ll200ll ; 
lHanll20o"5h . They are barely affected by each other unless the 
projected positions of the planets are close. On the contrary, 
in the central region of the magnification map the anomaly 
pattern induced by one planet can be significantly affected 
by the existence of another planet. 



3 DEVIATIONS DUE TO SECONDARY 
PLANET 

In Fig. 2, we show the fractional deviation of the magnifi- 
cation map obtained by fitting the binary lens model to the 
magnification map induced by triple lens systems, which is 
defined by 



Atri - Aj 

A hin 



(6) 



where A tr i and represent the magnification map gen- 
erated with the primary plus secondary planets and that 
obtained by fitting the binary lens model in which the 
secondary planet is absent, respectively. The inverse ray- 
shooting technique i s used to obtain the mag n ification map 
dSchneider fc Weiss! 1 19861 ; iKavser et all Il986j ; IWambsganssl 
Il997 ). Deviation maps are presented in terms of the lens 
position of the secondary planet both in angles between po- 
sition vectors of two planets 8 and in separations S2 (see 
Fig. 1). Each map of the fractional deviation is calculated 
as a function of the source position (£, rj). The coordinates 
are set so that the lens star is at the center. We concen- 
trate only on regions of j£j ^ 0.01 and |?7| S 0.01 for 
which high-magnification events are relevant. In generat- 
ing the magnification map induced by triple lens systems 
the primary planet is fixed to be located on the £ axis, 
(six, siy) = (—1.3,0.0). We set the mass ratio of the pri- 
mary planet to the lens star qi — 3 x 10~ 3 and that of the 
secondary planet 52 = 1 X 10~ 3 , corresponding to a Jupiter- 
mass and Saturn-mass planets orbiting a O.3A/0 star, respec- 
tively. Unless otherwise stated, in all computations through- 
out this paper we adopt typical values of distances for Galac- 
tic bulge events: D a \ = 6 kpc, D os = 8 kpc. We also assume 
the source radius 9+ to be 0+/0E = 0.001, which corresponds 
to a typical main-sequence star. Note that in the fitting pro- 
cedure parameters of the binary lens model composed of the 
a lens star and the primary planet are constrained to find 
their best fit values. To minimize the difference between the 
epsilon we generated using the triple lens and one we mod- 
eled by the binary lens model, the minimization process has 
been performed by the least square method with varying 
the mass ratio of the primary planet, q\, and its position, 
8\. Contours are drawn at the levels of e = ±1%, ±5%, 
±10%, and ±20%, and the regions of negative e are shown 
in blue and positive e in red. The color changes into light 
shade as |e| level decreases. 

As one may expect, distorted regions due to the sec- 
ondary planet are confined to the central caustic of the pri- 
mary planet due to the nonlinear interference between per- 
turbations produced by two planets. On the other hand, it 
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should be noted that when 8 = 0° there may exist a de- 
generate case where two-planet geometries will give rise to 
exactly the same magnification map with a sin gle planet of 
a sligh tly larger mass as previously noted by iGaudi et al.l 
(1993). In 

our particular example it can be seen in the case 
of S2 = 1.26 since the primary planet is located at (-1.3, 
0.0). We also note that the degeneracy between S2 and s^ 1 
can be seen due to the additional planet. Regions deformed 
by the secondary planet revolve as the angle between two 
planets varies, and become broader as S2 approaches unity. 
To explore the effect of q2, we show deviation maps for sev- 
eral S2 and Q2 in Fig. 3. We set the mass ratio of the primary 
planet to the lens star q\ = 3 x 10~ 3 and the angle between 
position vectors of two planets 9 = 30°. As in Fig. 2, con- 
tours are drawn at the levels of e — ±1%, ±5%, ±10%, and 
±20%, and the regions of negative e are shown in blue and 
positive e in red. From the maps, one finds that the lower 
limit of the mass of the secondary planet that may cause 
sufficient deformations in the magnification map is between 
1 x 10~ 5 < q2 < 1 x 10~ 4 , in which q2 — 1 x 10~ 4 corresponds 
to a IOMq planet orbiting a O.3M0 star. In other words, an 
Earth-mass planet can only be possibly detectable in a nar- 
row region of the parameter space unless the photometric 
accuracy allows one to examine the observational data in 
the levels of |e| ~ 1%. 



4 COMPARISON OF DETECTION 
PROBABILITY 

To 'discover' a planet, one may demand that A\ 2 between 
the light curve calculated by the planetary lens model and 
the observed data should be smaller than a carefully chosen 
critical value. Or, as another way to disclose a planet the 
Ax 2 -based criterion can be translated into so-called Gould 
& Loeb criterion, assuming a predetermined pho tometric 
accuracy and non- white noise (jGould fc Loebll 19921 ) . Gould 
& Loeb criterion considers deviations as a planetary sig- 
nal when a few observational points are consecutively de- 
viated from the single-lens light curve. In the current exer- 
cise, we are going to follow Gould & Loeb criterion. Thus, 
we will count a signal of planetary detection as successful 
if at least one point in the deviation map has the devi- 
ation amplitude larger than a certain threshold. We have 
defined the detection probability of the low-mass planet as 
follows: Firstly, for given q2, S2 and 9, we calculate the area 
in which the value of |e| in the deviation map, e.g., shown 
in Fig. 2 and 3, is greater than 5%. Bearing in mind that 
monitoring high-magnification events should be an efficient 
strategy of detecting 'super-Earth' planets in view of prac- 
tical purposes, we only consider this quantity in the range 
of |wo| ^ 0.01 as commonly employed. Having done that, 
we average the obtained area over the position angle in the 
range of 0° ^ 9 S 180°. Then, we normalize the averaged 
area with the area of |«o| ^ 0.01. In this way, we obtain the 
detection probability as a function of q2 and S2- Note that 
we repeat the same calculation for the binary lens with a 
planet that has the same mass as the secondary planet in 
the triple lens system, except that we construct fractional 
deviation maps given as 

4, • — A • , 

n bm ^single /—n 

£= A~- ', ' ^' 

^single 



where Abin and Single represent the magnification map gen- 
erated with the binary lens system having a planet and that 
obtained by fitting of the single lens model, respectively. 

In Fig. 4, we compare the probability of detecting the 
low-mass planet in triple systems with that in binary sys- 
tems, provided that the detection threshold is |e| > 5%. The 
detection probability of the secondary low-mass planet in the 
triple lens system is presented by grey-scale such that darker 
shade represents higher probability as shown in the scale- 
bar. We note that for a given §2 the probability in detect- 
ing the secondary planet becomes higher as S2 approaches 
unity, as expected in plots such as Fig. 2. We find that in the 
triple lens system the detection probabilities of the low-mass 
planet are ~ 50%, ~ 10%, and ~ 1% for q 2 = 10" 3 , 10" 4 , 
and 10~ 5 , respectively, if we only consider the secondary 
planet residing in the lensing zone, 0.6 < S2 % 1.6. For com- 
parison, we also present as dotted contours the probability 
of detecting the low-mass planet in binary systems. The de- 
tection probabilities for low-mass planets have been calcu- 
lated (Gould & Loeb 1992; Bennett & Rhie 1996; Griest & 
Safizadeh 1998; Gaudi, Naber, & Sackett 1998; Rhie et al. 
2000; Bennett & Rhie 2002; Rattenbury et al. 2002; Kubas 
et al. 2008). Interestingly enough, for a given 52 and S2 the 
detection probability of the low-mass planet in the triple 
lens system is very similar to that in the binary lens system. 
It makes sense when recalling that the magnification pat- 
tern due to the triple-lens systems is well approximated by 
the superposition of the magnifications due to the individual 
planets (Han et al. 2001; Han 2005). The fundamental rea- 
son why the detection probability of low-mass planets in a 
triple-lens system with a Jovian mass planet could be similar 
to the detection probability in a binary-lens system is that in 
most cases the binary superposition works well even in the 
central region. If the interference between caustics due to 
two planets destroys its original shapes then there is no rea- 
son that we end up with residuals implying a second planet. 
Having said so, another necessary condition that the detec- 
tion probability of low-mass planets in a triple-lens system 
with a Jovian mass planet could be similar to the detection 
probability in a binary-lens system is that the minimiza- 
tion should finds a solution for the Jovian planet accurately 
enough. 

In Fig. 5, we present the detection probability for the 
lower threshold of [g| > 1%, provided that photometric un- 
certainties of ^ 1% are achievable in the future. It is natural 
to find that the probability in the traditional lensing zone 
becomes higher. That is, we find that in the triple lens sys- 
tem the detection probabilities of the low-mass planet are 
~ 80%, ~ 40%, and ~ 5% for q 2 = 10" 3 , 10" 4 , and 10" 5 , 
respectively, if we only consider the secondary planet resid- 
ing in the lensing zone. It is also found that for high-mass 
planet of 52 > 10~ 4 corresponding to a WMm planet orbit- 
ing a 0.3Mq star there exist reasonably high probabilities 
in the broader range than the lensing zone, i.e., S2 < 0.6 
or S2 > 1.6. The detection probability of the low-mass sec- 
ondary planet in the triple lens system is also comparable to 
that in the binary lens system, except when q2 < 3 x 10 -6 
corresponding to a O.3M0 planet orbiting a 0.3Mq star. In 
that parameter space the detection probability in the binary 
lens system is somewhat higher than that in the triple lens 
system. 
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5 CONCLUSION 

We have computed the detection probability of the low-mass 
planet specifically for high-magnification events of the triple 
lens system, motivated by the fact that the central caustic 
is distorted by any companion to the lens star. Having done 
that, we have compared the detection probability of the low- 
mass secondary planet for high-magnification events of the 
triple lens system with that of the same low-mass planet 
but in the binary lens system. It should be stressed that 
the detection probability of the low-mass planet for high- 
magnification events of the triple lens system given in this 
paper is a kind of relative probability in the sense that our 
criteria on the discovery of a planet is actually closer to one 
of signals of the planet. The detection of signals of a low- 
mass planet may not be sufficient for the discovery of such a 
planet. Having compared results of two cases with the same 
criteria, however, probabilities given here are still worthy in 
lending emphasis to the light curves of the triple lens system. 
Our main findings are as follows: 

(1) In the triple lens case where the secondary planet 
resides in the lensing zone, the detection probabilities of 
the low-mass planet for the high magnification events are 
~ 50%, ~ 10%, and ~ 1% for q 2 = 10~ 3 ,10~ 4 , and 10~ 5 , 
respectively, when the detection criterion is jej > 5% where 
e is the deviation in magnification. 

(2) When the detection criterion is |e| > 1%, those 
probabilities increase to ~ 80%, ~ 40%, and ~ 5% for 
q2 = 10~ 3 ,10~ 4 , and 10 -5 , respectively. For high-mass 
planet of q-z > 10 -4 , there exist reasonably high probabilities 
outside the usual lensing zone, S2 < 0.6 or S2 > 1.6. 

(3) For a given qi and S2, the detection probability of 
the low-mass planet in the triple lens system is compara- 
ble to that in the binary lens system. Therefore, it is quite 
necessary to develop an efficient algorithm search for 'super- 
Earth' planets in the paradigm of the triple lens model as 
well as of the binary lens model. 
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Figure 1. Geometry of triple lens system composed of a lens star, a primary planet and a secondary planet. In the upper panel, the 
coordinate (£, r/) are centered at the position of the lens star, and 6 is the angle between position vectors of two planets (qi, 52). All 
lengths are normalized by the radius of the Einstein ring. The gray scale represents the magnification and brighter tone represents 
higher magnification. The lower panel shows the zoom of the boxed region of the high magnification within uq ^ 0.01. 
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Figure 2. Fractional deviation in magnification of triple lens system from the binary lens system. Deviation maps are presented in 
terms of the lens position of the secondary planet both in angles 9 between position vectors of two planets and in separations S2, 
each of which is calculated as a function of the source position (£, if). The separation of the low-mass planet and the angle between 
two planets are marked on the right side of each row and the top of each column, respectively. The coordinates are set so that the 
lens star is at the center. The primary planet is fixed to be located on the £ axis, (si x , si y ) = (—1.3, 0.0). We set the mass ratio of 
the primary planet to the lens star qi = 3 X 10~ 3 and that of the secondary planet q2 = 1 X 10~ 3 . Contours are drawn at the levels 
of e = ±1%, ±5%, ±10%, and ±20%, and the regions of negative e are shown in blue and positive e in red. The color changes into 
light shade as |e| level decreases. 
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Figure 3. Similar maps as Fig. 2, except that we set the angle between position vectors of two planets 6 = 30° and vary the mass 
of the secondary planet. As in Fig. 2, contours are drawn at the levels of e = ±1%, ±5%, ±10%, and ±20%, and the regions of 
negative e are shown in blue and positive e in red. 
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Figure 4. Probabilities of detecting the low-mass planet in binary systems and in triple systems. The detection probability of the 
secondary low-mass planet in the triple lens system is represented by grey-scale, and is drawn such that darker shade represents 
higher probability as indicated in the gray index. For comparison, we also present as dotted contours the detection probability of 
detecting the same low-mass planet if it is in a binary systems. Probabilities are calculated such that the value of |e| in the deviation 
map is greater than 5%, considering only \uq\ ^ 0.01 events. 




Figure 5. Similar maps as Fig. 4, except that probabilities are calculated such that the value of |e| in the deviation map is greater 
than 1%, considering only |uq| < 0.01 cases. 



